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@ Reduction of image artifacts from support structures In tomographic imaging. 



® An apparatus for reducing Image artifacts 
caused by over-ranging or clipping of the data 
collected in a tomographic scan using a patient 
support models the density profile of the pa- 
tient support with a polynomial. Coefficients of 
the polynomial are stored m a look-up table 
according to gantry angle and used to estimate 
the over-range data. The estimated data is blen- 
ded with the in-range data by convolving an 
over-range mask indicating which data is in- 
range and which data is over-range with a box 
car convolution kernel which produces a 
trapezoidal weighting mask. This weighting 
mask is multiplied by the estimated data and 
then summed to the in-range data to correct for 
the over-range portions of the data coBected. 
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Background of the invention 

This invention relates to tomographic Imaging systems, such as x-ray computed tomography. More spe- 
cifically, the Invention relates to a method for reducing Image artifacts caused by structures used to support 
an imaged object 

Tomographic Imaging systems, as considered herein, are imaging systems which produce a tomographic 
or "slice 0 picture of an object such as the human body. Such tomographic systems collect a sertes of projec- 
tions 0 at various angles around the body, each projection made up of measurements of radiation emitted from 
the body. The radiation may be that transmitted through the body by an external radiation source, such as an 
x-ray tube, as in x-ray computed tomography ( n CT) or the radiation may be that emitted from the body from 
internally placed radk>-pharmaceutlcaJ Isotopes, e.g. T c ®> used In Single Photon Emission Computed Tomog- 
raphy (SPECT). 

In an x-ray CT machine, an x-ray source is collrmated to form a planar x-ray beam within an x-y plane of 
a Cartesian coordinate system. The beam is transmitted through the imaged body and received by a generally 
linear detector array also within this x-y plan®. 

The detector array is comprised of a plurality of adjacent detector elements each receiving radiation along 
a single W from the focal spot of the x-ray source to that detector element Each detector element produces 
a signal Indicating the attenuation of the x-ray beam along that ray by the imaged body. The detector elements, 
the signal from each detector element, and the ray of the fan beam are all generally referred to as a °channei° 
of the projection, any ambiguity generally being resolved from the context of the use. 

In one common embodiment, the planar x-ray beam Is a "fan beam" radiating from a focai spot and the 
detector array has its elements organized in an arc of constant radius about ths focal spot The x-ray source 
and detector array may be rotated on a gantry around the imaged object so that the fan beam intersects the 
imaged object at different angles. 

A number of projections are acquired at different gantry angles to form a tomographic projection set The 
acquired projection set is typically stored in numerical form and may be Reconstructed" by mathematical tech- 
niques to yield a slice image. The reconstructed images may be display on a conventional CRT tube or may 
be converted to a film record by means of a computer controlled camera. 

The natural evolution of x-ray CT has led to the development of higher powered x-ray sources. Such x- 
ray sources produce increased x-ray flux which is desirable for two reasons: First, increased flux improves the 
signal-to-noise ratio in the resulting image, for example, by minimizing the effect of noise. Second, Increased 
flux parmits the acquisition of the projection set during a shortened scanning time. Decreasing the scanning 
time improves patient comfort and helps minimize image artifacts caused by patient motion. 

The use of increased x-ray power In an x-ray CT system creates the potential of overwhelming the system's 
detector signal chain. In particular, the analog-to-digital converter f ADC 11 ) associated with the CT system's 
data acquisition system ("DAS 0 ) may be driven over its range. Such an over-range condition artificially limits 
the signal from any over-range channel to the maximum value of the ADC and causes the data of these chan- 
nels to be effectively lost 

For the case of a patient scanned by a CT system, the over-range channels wfll typically be those channels 
which receive the peripheral rays of the fan beam, e.g. those near the outer edges of the patient and beyond 
These areas may include support structures such as the patient table or a headhoider end the loss of the chan- 
nel data for these structures. Although it is relatively unimportant to produce accurate images of the support 
structure, Incomplete data in this area will cause artifacts throughout the entire reconstructed image. This 
spreading of the effects of locally erroneous projection data is caused by the frequency domein filtering implicit 
in the fmage reconstruction process. 

The problem off over-range channels is s Imflar to t he problem created by the truncation of projections when 
the imaged object extends outskte of the fan beam and detector range. Such truncation may occur, for example, 
when the patient's arms are outside of the fan beam for some projections and in the fan beam for other on> 
factions, 

U.S. Patent 4,550,371, incorporated by reference and assigned to the same assignee as the present ap- 
plication, provides a means for correcting a truncated projection by calculating the moments of ail projections. 
Unfortunately, this method relies on the assumption that the majority of the projections are not truncated, an 
assumption rich does not hold, in general, for projections having over-range channels. 

One proposed method of correcting for the effect of over-range channels involves independently meas- 
uring the thickness of the Imaged object along the rays associated with the over-range channels and substi- 
tuting an attenuation value for those rays based on assumption of constant density of the imaged object along 
those rays. Although some success has been achieved with this approach, the requirement that the thickness 
of the body be measured independently is cumbersome and commercially impractical. Further, such meas- 



urements, based on simple models of the body being imaged do not take into account the complex attenuations 
caused by proximate patient supporting structure. 

Ideally, any system for correcting for erroneous projection data caused by over-ranging should operate 
quicCdy enough so as not to significantly delay the production of tomographic images from the projection set 
ff Preferably, the correction system should allow correction to begin as the projections are being acquired and 
should be susceptible to parallel processing in an array processor or the like, such processors as aire commonly 
used for the reconstruction of tomographic images. 

Summary of the Invention 

10 

In one aspect, the invention provides a compensation circuit for use In a tomographic scanner, the scanner 
employing a fixture adjacent to a body to be scanned and measuring a set of intensity values each related to 
the Intensity of received radiation along one of a set of adjacent rays through both the scanned body and the 
fixture from a given gantry angle, the intensity values indicating the intensity of the received radiation up to a 

16 maximum, and not indicating the intensity of the received radiation above that maximum, the sat of intensity 
values being combined with other sets of intensity values and reconstructed to product an image, the circuit 
comprising a threshold detector receiving the set of intensity values for Identifying the set of intensity values 
as in-range intensity values corresponding tp received radiation not above the maximum, and over-range In- 
tensity values corresponding to received radiation above the maximum; an model generator receiving t he value 

20 of the gantry angle for producing estimated over-range intensity values corresponding to the intensity of radi- 
ation along rays associated with those values through the fixture alone at the gantry angle; and a combiner 
for combining the in- rang a intensity values with the estimated over-range intensity values prior to reconstruc- 
tion of the intensity values into an image. 

In a further aspect, the invention provides In a tomographic scanner measuring a set of intensity values 

as each related to the intensity of received radiation along one of a set of adjacent rays through both an imaged 
object and a fixture , from a given gantry angle, the fixture being adjacent to the body to be scanned, tho in- 
tensity values indicating the intensity of the received radiation up to a maximum, and not indicating the intensity 
of the received radiation above that maximum, the set of intensity values being combined with other seta of 
intensity values and reconstructed to produce an image, a method of correcting for over-range comprising gen- 

30 erating an over-range mask identifying the set of intensity values as in-range intensity values corresponding 
to received radiation not above the maximum, and over-range intensity values corresponding to received ra- 
diation above the maximum; determining estimated over-range Intensity values corresponding to the Intensity 
of radiation along rays associated with those values through the fixture alone at the gantry angel; and com- 
bining the in-range intensity values with the estimated over-range intensity values prior to reconstruction of 

ss the intensity values into an Image. 

The present Invention provides a means for correcting over-range data which results in the truncation of 
projections, specifically where the truncated data includes portions of a patient support structure or the like. 
The means requires as an input only the projection data from the projection being corrected and may be per- 
formed in parallel on the data processing hardware typically associated with tomographic systems. 

to Specifically, the channels of the projection data are received by a threshold detector which identifies them 
as in-range channels orover-range channels based on the maximum value associated with the data acquisition 
system. A model generator, programmed with the characteristics of the patient support, produces estimated 
over-range channels for the patient support Acombiner receives the estimated over-range channels end com- 
bines them with the in-range channels prior to the projection set being reconstructed into an imege. 

4$ It is one object of the Invention to provide an accurate estimation of the values of the over-range channels 
In a projection set without the need to scan the data of the other projections of that set. The model of the patient 
support uses gantry angle to deduce the contribution of the patient support to the present projection. 

The model of the patient support may employ a function generator which receives a set of coefficients 
baaed on gantry angle to model the patient support by means of a polynomial. The polynomial may employ a 

so pieca-wise Gaussian basis function. 

Thus, It is another object of the invention to provide an accurate model of the projections of a complex 
patient support without the need to store extremely large amounts of data, such as that which would be required 
to store the actual projections of the patient support. The use of a piece-wise Gaussian function as the basis 
function for the polynomial allows the patient support to be adequately characterized with as few as six coaf- 

66 fici© nts per gantry angle. 

The combiner may comprise a mask generator for generating a binary over-range mask Indicating which 
channels are over-range and which channels are in-range. A box car generator may be used to produce a con- 
volution kernel which may be convolved with the mask generator signal to produce a convolved signal. The 
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estimated over-range channels may be multiplied by the convolved signal and summed with t he in-ranga chan- 
nels to effect the combination. The convolved signal may be the convolved over-range masCc and kernel dipped 
and scaled by a predetermined value to control the space between the estimated over-range channels and 
the In- rang® channels. 

5 Thus, It Is another object of the invention to provide a correction circuit suitable for Implementation In par- 
allel on an array processor or the like. The operations of convolution* multiplication, dipping and scaling are 
all readBy performed in parallel as required by such processors. 

The foregoing and other objects and advantages of the invention will appearfrom the following description. 
In the description, reference is made to the accompanying drawings which form a part hereof and In which 

10 there Is shown by way of illustration, a preferred embodiment of the Invention. Such embodiment does not nec- 
essarily represent the full scope of the Invention, however, and reference must be made therefore to the daims 
herein for interpreting the scope of the invention. 

Brief Description of the Drawings 

15 

Figa. 1(a) and (b) are front and side views, in elevation, of a CT gantry showing the relative positions of 
an x-ray source, detector array and fan beam about a patient's head held within a headholder; 
Fig. 2 is a block diagram of a CT control! system associated with the gantry of Fig. 1 and useful for practidng 
the present invention; 

Fig. 3 is a block diagram showing the processing of the data acquired by the CT system according to the 
present invention; 

Fig. 4(aHd) are graphical representations of signals producod In the present invention by an example ob- 
ject held in the headholder of Fig. 1; 

Fig. 5 Is a graph showing the attenuation versus channel for the hsadholder of Fig. 1 at a gantry angle of 
0 degrees; 

Fig. 6 is a figure slmHar to that of Fig. 5 showing the density prof ile of the headholder of Fig. 1 at 45 degrees; 
Fig* 7 is a graph of the basis function employed to model the projections shown in Rga. 2 and 3 being a 
piece-wise Gaussian function; 

Fig. 8 Is a perspective view of the headholder of Fig. t ; and 

Fig. 9(aH<J) are graphical representations of an over-range mask and signals generated therefrom for cor- 
recting the signals of Fig. 4(c) per the present Invention. 

Detailed Description of the Preferred Embodiment 

35 CT System Hardware 

Referring to Fig. 1 , a CT gantry 16, representative of that used with a "third generation* CT scanner, hdds 
an x-ray source 10 producing a fan beam of x-rays 24. The fsn beam 24 is directed through a patient 12, posi- 
tioned near a center 19 of the gantry 16, to be received by a detector array 18 also attached to the gantry 16. 
40 The patient's head is supported by a headhdder 1 5, the latter producing a small but significant attenuation of 
the x-rays of the fan beam 24 . 

The gantry 16 rotates within an x-y plane of a Cartesian coordinate system, termed the imaging plane, 
which is generally the same plane as that of the fan beam 24. 

The detector array 18 is comprised of a number of detector elements or "channels 0 26 positioned adjacent 
45 to each other within the imaging plane to subtend the fan beam 24. The channels 26 receive and detect radi- 
ation passing from the x-ray source 10, to produce a plurality of channel signals each associated with a par- 
ticular channel 26. At a given orientation of gantry 16 about patient 12, signals for approximately 600 channels 
may be acquired, representing a delated picture of the line integral of the attenuetion of the fan beam 24 by 
the patient 12 at that angle. A gantry angle of zero is defined as that angle where a prlndple ray 28, centered 
so in the fen beam 24, is directed vertically downward from the x-rays source 1 0. 

The x-rays of the fan beam 24, Immediately after leaving x-rays source 10 and prior to being received by 
the detector array 1 8, air© filtered by a spectral filter 30 which filters out the lower energy x-rays from the fan 
beam 24. The fan beam 24 then passes through a bow tie filter 32 having a profile that produces an attenuation 
in the fan beam 24 complementing that which would be produced by a cylinder of water placed at the center 
55 ig of the gantry 16. The purpose off the bow tie filter 32 is to reduce the range of Intensity values received by 
the detector channels 26 for a typical patient 12 and hence to aOow for an increase in sensitivity of the detector 
array 18 and its associated circuitry. This increased sensitivity also, however, increases the chance of a chan- 
nel being over-range when the patient 12 differs significantly from the water cylinder model. 
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Th© bow tl© filter 32 Is followed by an aperture 34 which forms fan basm 24 and may bo used to correct 
the position of the fan beam 24 with respect to the surface of the detector array 1 S as described generally in 
U.S. Patent 5.054 ,041 issued to the same assignee as that of the present application and incorporated heroin 
by referanc®. 

For a given patient 12, the channels 28 may be roughly divided into three groups: reference, over-range, 
and in-range. Reference channels 20 of the detector array 18 are those intended not to be occluded by the 
patient 12 or headhoSder 15 and may serve the function of calibrating the projection data for variations in the 
x-ray flux from x-ray source 10, and sent® further to permit automatic alignment of the fen beam 24 on the 
detector array 18. Over-range channels 21 of the detector array 18 are those channels within e given projection 
which, although possibly occluded by the Imaged object 12, generally receive x-rays having so little attenuation 
that the ADC. used to digitize the signals of t hose channels, is over-ranged. This over-ranging wfll be described 
further below. And finally, In-range channels 22 of the detector array 18, frequently but not necessarily near 
the center of the fan beam 24, are those In a given projection which are sufficiently attenuated by the Imaged 
object 12 so as not to over-range the ADC used to digitize the signals from these channels. 

fin general, soma ovor-rang© channels 21 will be present for substantially all of the projections ta&en of a 
convex patient 12 subtending less than the entire fan beam. 

Referring to Rg. 2, control circuitry for a CT imaging system suitable for use with th© present invention 
includes a number of functional blocks 48. A data acquisition system 62 such as is generally understood In the 
art, is connected to the detector array 18 and comprises a sampling means (not shown) for sampling the signals 
from each of the channels 26. An analog to digital converter {"ADC*) (not shown) converts the sampled analog 
signals from each sampled channel 26 to a digital value for processing by later circuitry. The ADC has a finite 
range and generally a trade-off must be established batwean the range of th© intensity signals from each chan- 
nel 26 which may be correctly digitized, and the resolution of the digitization process. Both considerations 
are important and a reasoned trade-off will necessarily allow certain situations where the rang® of the ADC 
wiil be exceeded, in the case of an over-range, the ADC will simply output its maximum value, regardless of 
how much greater th© signal from the channel 26 is than this maximum value, 

Aradiotranslucent table 53 supports the patient 12 and the headholder 15, the latter which Is typically fixed 
to the table 53. The table 53 may be moved t hrough the image plane to align the slice of interest of the patient 
with the image plane, and may be raised or lowered to center the patient 12 within the opening of the gantry 
16. The movement of the table is accomplished by motors (not shown) controlled by table motor control 51. 
The table motor control 51 also generates a value H Indicating the height of table 53 with respoct to the too- 
center 19. 

An x-ray control 54 provides power and timing signals to the x-ray source 10 with regard to the position 
of gantry 16 to acquire the projections. Gantry motor controller 56 controls the rotational speed and position 
of the gantry 16 and provides gantry angle Information 0 to the DAS 62 and the x-ray control 54 to permit 
accurate timing of the projections. 

The image reconstruct 68 is a special purpose computer, such as an array processor. capaWe of very 
rapid parallel processing or "pipelining* as is necessary to produce images from the large amount of projection 
data. Array processors suitable for use as the image reconstructed 68 are commercially available from a variety 
of sources. The image reconstructs 68 receives the sampled and digitized signals from the channels 28 of 
the detector array 1 8 via the DAS 62 to perform high speed image reconstruction according to methods known 
in the art. 

A computer 60 coordinates the operation of the DAS 62, the table motor control, the x-ray control 54, and 
the gantry motor control 56 and works in conjunction with image reconstructor 68 to reconstruct tomographic 
images from the set of projections acquired by the scanning process. The computer 60 receives commands 
and scanning parameters via operator console 64 rich fe generally a CRT display and keyboard which allows 
the operator to enter parameters for the scan and to display the reconstructed slice images and other infor- 
mation from the computer 60. A mass storage device 66 probldes a means for storing operating programs for 
the CT imaging system, as well as image data for future reference by the operator. 

Each of the above elements is connected to Its associated elements on the gantry 16 via slip rings 50 to 
permit continuous rotation of the gantry 16. 

Operation of the Correction Process 

Referring now also to Figs. 4(a) and (b), for a given position of gantry 1 6. the DAS 62 will produce a raw 
projection 36 comprised of t he stg nals from each of t he channels 26 of t he detector array 1 8. This raw projection 
36 has an in-range portion 22 1 corresponding to the in-range channels 22 and an over-range portion 2V cor- 
responding to the over-range channels 21. For certain gantry angles, the over-range portion 21 may Include 
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some in-range channels associated with t ha headholder 1 5 which may provide sufficient attenuation to prevent 
over-ranging. The channel values of the over-range portion 2V are limited or clipped to the ADC maximum 
value 39, as has been described above, to produce clipped data 38. 

Referring momentarily to Fig. 4(a), the in-range portion 22\ In this example, corresponds to x-rays passing 

5 generally through the midsection of both the patient 1 2, and the bow tie filter 32 whereas the over-range por- 
tion 2T corresponds to x-rays intersecting only the edge and hence relatively little thickness of the patient 12 
and the headholder 1 5 and th© edga of the bow tie filter 32. The increasing attenuation of the bow tie filter 32 
of the x-rays of the over-range portion 21\ ultimately brings the raw projection signal 36 bade In-range for chan- 
nels at the extreme edge of the detector array 18. However, at the point where the x-rays of the fan beam 24 

10 first pass on either side of the patient 12 without attenuation by the patient 12, the thickness of the bow tie 
filter 32 and the headholder 15 will often be insufficient to prevent clipping of the raw projection signal 36 at 
over-range portion 21'. 

It will be understood that In situations other than that illustrated, the over-range pardon 21 ' may be toward 
the mid-section of the detector array 18, for example if the Imaged object 12 is substantially displaced about 
is the isocenter 1 9. Thus t he In-range portion 22* and over-range portion 21 ' may generally vary in their positions 
with respect to each other and the detector array 18. 

Referring now to Figs. 4(b) and 4(d), the clipped data 38 of the raw projection 38 Is readily Identff bd as 
those portions of the raw projection 36 equal to the ADC maximum 39. Accordingly, an over-range mask 40 
may be simply generated by comparing the raw projection 36 to the ADC maximum 39 and setting the over- 
range mask to °zero° for those points of the raw projection signal 36 substantially equal to the ADC maximum 
39 and setting the over-range mask value to D one q for ail other points. 

Reforrlng now to Fig. 3, this thresholding operation la represented by process block 42 which may be im- 
plemented by the array processor of the image reconstructed 68. 

The raw projection 36 of Fig. 4(b) is next adjusted for the effects of the bow tie filter 32 and channel to 
channel verlation In detector array 18 at an air correction step, indicated by process block 82, to produce cor- 
rected data 92. The correction of process Wock 82 is simply a point by point multiplication of the raw projection 
data 36 times a function, representing the compliment of the attenuation of the bow tie filter 32 and the channel 
to channel variation in the detector array 18, to effectively remove the contribution of the bow-tie filter 32 from 
the corrected data 92. This correction is followed by other calibration steps such as beam hardening correction 
» etc. 

Referring to Fig. 4(c), the corrected data Is next log adjusted by taking the negative of Its logarithm to pro- 
vide a density profile 44. Density In this case refers not to mass per unit volume but simply to the total amount 
of attenuating material of the imaged object 12 along the x-ray beam associated with a particular channel 26. 
The logarithmic correction Is required because the attenuation of x-rays (or any radiation) by a medium Is oen- 
erally exponential, as illustrated by the following formula: 

I - to® - \v*m 

where lo is the x-ray intensity before it passes through the object having attenuation coefficients u given over 

the x-y plane as n(x,y), I is a distance along the x-ray path length and I is the x-ray intensity after the object 
The negating operation simply reflects the fact that decreasing Intensity represents increasing density. 
40 This correction process is represented in Fig. 3 by process block 47. 

Referring still to Fig. 4(c), the density profile 44 includes erroneous data In the areas 46, associated with 

the raw projection data of over-range portion 21' of Fig. 4(b). The correct data in these areas 46 are shown 

by the dotted lines 48 flanking the central undipped data of the density profile 44, the latter undipped data 

associated with the raw projection data of in-range portion 22' of Fig. 4(b). 
45 Although the dipped channels 46 of the density profile 44 are at the edge of the density prof Be 44. they 

win create artifacts throughout the reconstructed image as a result of a convolution of the density prof ile 44 

in the image reconstruction process as will be described. 

Convolution has the effect of spreading the contribution of each channel of data In the density prof U 44, over 
the entire reconstructed image. Accordingly, it is necessary to correct even peripheral areas of clipped chan- 

50 nets 46 prior to image reconstruction. 

The present invention provides corrected data for dipped channels 46 by considering two factors. The first 
is the contribution of the patient 12 to the data of these dipped channels 46. The second is the contribution 
of the headholder 16. The present invention provides estimates of bothof these contributions to generate sub- 
stitute values for the data in the dipped channels 46 without reference to previous or later projections and with- 

55 out external measurements of the patient 12. 
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The Patient Contribution 

The estimate of the contribution of the patient 12 is accomplished by employing a simplified model of the 
imaged object and fitting that model to the In-range channels for that projection. This process Is the subject 
6 of co-pending application serial number (15CT03712) filed on even date herewith and entitled: "Over-range 
Image Artifact Reduction in Tomographic Imaging" as now described. 

The use of a model to calculate the patient contribution is based on the recognition that the slope of the 
data of the projections at the start of the areas of clipped data 46 will be directly proportional to the size of 
the object 12. This relationship may be illustrated in the following description where ft will be assumed that 
10 the rays of the fan beam 24 are parallel rather than fan beam shaped for the purpose of clarity. The extension 
of this description to the case of e fan beam is a straightforward geometric transformation that wiP be under- 
stood to those of ordinary skill in the art It will also be assumed that the imaged object 12 is a simplified geo- 
metric solid, preferably as a cylinder having radius R and a uniform attenuation coef ffcient p. This requirement 
will be relaxed later. 

is The density profile of the cylinder is described by the following equation: 
P<x) - jjjWR 2 - X 2 if X S R 



20 



28 



30 



50 



[0 if X > R <*> 



The term 2Vr 2 - x* is simply the path lengt h along a chord through a cylinder of radius Rasa function 
of its x-coordinate, the x-coordinata being measured along the detector array 18 with a value of x » 0 at the 
center of the detector array 18. For the parallel case, the detector array 18 will be flat rather than curved. 

For this cylinder, an over-range of a channel In the density profile 44 will occur for channels where the 
path length p(x) Is less than a predetermined value C. The value C depends on the attenuation coefficient u 
of the material of the cylinder, the strength of the x-ray source 10 after passing through the filter 30 and bow 
tie 32. and on the ADC maximum value 3$. Taking the derivative of equation (1), and evaluating It at C provides 
the slope 58 where over-range starts: 

dpftO „ VR*-C* „ M 

ir ■ — c — atx " c (2) 

Since C will be much less than R for most clinical situations, equation (2) can be simplified to- 

Equation (3) indicates a linear relationship exists between the radius R of the imaged object 12 and the 
slope 58 of the density profile 44 at the starting point 52 of the clipped channels 46 (all shown In Fig 4(c)) 
where over-range starts. The fact that this relationship is linear in this example greatly simplifies the estimation 
of the data of the clipped channels 48, it will be understood that other geometries of the imaged object 12 may 
be substituted for the cylindrical geometry If increased complexity of calculation performed above may be tol- 
erated. Important, primarily, is that equation (3) suggests that slope 58 will be a function of R for a range of 
Imaged objects 12, and as will be explained below, R may be deduced from the data of a single projection. 

In order to take advantage of this functional dependence of the slope 58 of the data In region 46. on the 
value of the radius R. one must have an estimate of R. This estimate of R Is obtained as follows: If. for the 
45 assumed parallel geometry, the attenuation of each volume element of the Imaged object 1 2 in the imaging 
plane is represented by M (x,y), a projection at a given gantry angle 9 may be represented by the following equa- 



p (*/t)« / / /t(x,y)3(xco«* ♦ yainf-tjdxdy (4) 



55 where S is the delta function and simply provides a notationaily convenient way of generating the line integrals 
through n(x,y) of the projection at gantry angle d, and t Is a function of x and y related to the distance of each 
ray of the projection from the center ig. along a line normal to the path of the ray. The value t corresponds 
roughly to the position of the channels 26 within the detector array 18. 
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If we Integrate this equation (4) with respect to the variable t, as would be done in practice by summing 
each value of the density profile 44 together, we obtain; 



5 M " / / / ^( x 'y)*( xc ° 8 ** y«in5-t)dxdydt 

• 09. CO -oo 
OO 00 

» - J J M<x#Y)dxdy (5) 

-09 -00 

Equation (5) is no longer a function of 0 indicating that the sum of all the data in the density profile 44 
is depends solely on the attenuation coefficients u.of the object and is essentially Independent of the gantry angle 
8 and thus of the particular projection acquired. Although this Is not true for fan beam geometry, It has been 
determined thet the variations in this integral for a fan beam case are small and thus the integral of equation 
(5) can be considered constant from projection to projection even for fan beams. 

The total attenuation given in equation (5) may be used to estimate the value R. For a cylindrical imaged 
20 object the relationship is according to the following formula: 

(6) 

In the preferred embodiment, however, it is assumed only that the slope 58 is functionally related to the 
value of M and the functionality is determined empirically. Experiments have shown that the relationship be- 
25 tween slope 58 and M is roughly linear in the regions of interest and therefore a simple constant of proportion- 
ality suffices. This constant of proportionality is obtained from measurements of actual patients. In practice 
different constants of proportionality are used.for Images of a patient's head as opposed to images of the pa- 
tient's torso. 

Referring to Fig. 3, this summation to determine M Is shown by process block 74. The slope calculation 
^ of equation (3) is represented by process block 72 of Fig. 3. 

Referring again to Fig. 4(c), the data of the density profile 44 at areas 46 may be conformed to the value 
of the slope 58 produced by the slope calculator 72 simply by searching the density profile 44 for the starting 
points 52 of the clipped channels 46. Working forward and backward from these starting points 52, the density 
of the starting points 52 may be progressively decreased according to the calculated slope 58 to generate new 
density values for the clipped channels 46. 

Preferably, however, this correction process employs a method compatible with the pipelining capabilities 
of the array processor in the image reconstructor 68 as will be described below. 
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The Headholder Contribution 

Referring to Fig. 8, the headholder 15 comprises a U-shaped trough of Plexiglass 96 having a bottom foam 
pad 98 for cradling the patient's head. 

The headholder 15 provides an invariant contribution to the density profiles 44 of a tomographic scan em- 
ploying that particular headholder 15, in contrast with the variations expected In the particular patients 12. 
Thus, the headholder 15 may be precisely characterized without resorting to the approximations of a simplified 
geometric model as is done with the patient 12. 

Unfortunately, precise characterization of the headholder 15 is daunting. As shown in Figs. 5 and 6, the 
densfty profiles 44 of the headholder 15 at various gantry angles show considerable variation. Even for a single 
density profiles 44 of the headholder 15, it is evident that the density profile 44 may not be approximated with 
a simple geometric model as was done with the patient 12. For example, the profile density 44 at a gantry 
angle of 0 degrees shown in Fig. 5, includes two relatively sharp spikes flanking a plateau. In contrast, the 
denstty profile 44 of the headholder 1 5 at the gantry angle of 45 degrees is completely asymmetric. 

The discontinuous character of both density prof fles 44 of Figs. 5 and 6 suggests that accurately digitizing 
and sampling these density profiles 44 would require a large amount of stored data. Further, the variation in 
t he density prof Pe 44 of t he head holder 1 5 wit h gantry angle requires that multiple density profiles 44 be stored , 
one for each gantry angle, if this density profile 44 data is to accessible for error correction in the same time 
frame as is the contribution of the patient 12 described above. 

In order to avoid the burdensome amounts of data needed to store the density profiles 44 of the headholder 
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15, the present invention models the headholder 15 with a polynomial. Only the coefficients of the polynomial 
are stored for each gantry angle, substantially reducing the amount of stored data. 

Given the general symmetry of the headhoJder. for example, as Illustrated in the density profile 44 of Flo 

IXS^riHSTt "* ^ IT.* * PO,ynomW ***** 8 ™P basis function <*i 

^/r, !'", ) WOuW bB ,nappr °P riat8 - ,nste ad customized basis functions were generated, one for diar- 

mXSS^T^ ° f Profi ' 8 44 "* ° therfor Ch8f8Cterfein0 th6 r, « ht haraJ ^tion 

44 rjS!! 7 ? basls 1 func,ton 100 for ^composing the left hand portion of the density profile 
TJl ™„ 1, 5 ' 8 ° Pfeca " wl8e comb,naUon of Gaussian functions, normalized to values between o 
and 1 and reaching an asymptote of approximately 0.35 aftarfiiatrisingto a peak value of 1. This basis function 

Z ^ .Tf °! th8 <ten8My Pr0,il8 44 ° f ,hB h88dh0,dar 15 as shown in Fig. 5 but is suS 

STlSiSSr 8 44 8h0W ° 10 Fi9t * Th8 Mcond 08818 <"°t shown) h «£* 

the left to right mirror image of the basis function shown in Fig. 7. 

The headholder for a given gantry angle may be now characterized by the following equation. 



<cr a <0)8k <■ ffgtflfb) (7) 



where 0 is the gantry angle and %) is the density profile 44 of the headholder at that gantry angle and 

£LHL ?& l afe ., th f ' 8ft ^ riflW Da8i8 fUnCttons re8 P 8Ctlve| y- The coefficients 0,(9) and JU8) may be 

^SmSST* by w mea6Urtnfl 4 nUm ° er * prtrfIles 44 * a headho,d ^ S without apatSli 
h ! coa ( ffK:le " ts to lha actua » Pylons entails dividing the density profiles 44 into a lefi 

S hS 2?J rf.iT 1 °V helr ,0ta ' ,9nflth and U8infl tha 088,8 100 of fJ. Vto charectenle the 
h«L m ?! , 8 b , a8i8 fUnC,i ° n lh8t 18 lhe m,rror ,ma 8 8 * *™*on 100 to charactered right 

It has been determined that each density profile 44 of the headholder 1 5 at a given nantrv anole can h« 

chants EASTS' pr0f,,e f 4 ^ tha "aadholder 15 Is next shifted with respect to the corresponding 

« a-ZJi^L ^ S 8nd hCnCe ,hfi P 08 ' 110 " of tha actual headholder 15 with respect to thTieocJnlar 

££^££2??™ *? f0re9 ° in8 de8Criptlon ' to th08e 8kfl » ,n tha ^ the anTnTof 

auch shifting wis depend also on the gantry angle 0 according to well known trigonometric formulas TNtaaMf? 
Ing of the generated density profile 44 of the headholder 15 aligns the mnMd^ah2^J^^ 
ectoal position of the headholder 15 for each protection as attach to the S 53 * ^ ^ 

The process of generating the headholder projection used in the present invention from the oantrv ancle 
fe shown a6 ^ b ock 102 of Fjfl 3 Proca88 Wock , Q3 £ y SK^SS* TOfc 

gerrtry angle 8 is used to look up in a lookup table (not shown) the necessary *%S^£!a2^ 
put into the equation 7 to generate the necessary density profile 44 coe ^^ wh.ch ™ then 



Combining the Patient and Headholder Contributto; 



ins 



it is necessary that the patient end headholder contributions to the data of areas as d««»«ih~t -k^ ^ 
combined smoothly with , he in-range channeis to avoid introduction"^ 

construct.on process. This may be accomplished by a weighting function used to weight the Sient end 
h*de« 

oen^ih-TS ' nV9nti0n ' the Wel9hl,n9 8nd COmbln,nfl proce88 18 ec^Pltehed so as to be amenable to 
computation on the array processor 68. According.*, referring to Fig. 3. the slope calculated by Ihe^ooe cat 

£2 a ST rfl^i 3 ^ re8PeC< 10 th8 »"«»*'" of the patient 12 described above, is ^ to S 
duce a boxcar kernel 75. The box car kernei 75 is a binary signal having a values of 1 for a pu W wTdlh^ 
chosen so that t = 

slops 
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The box car kernel 75 and the over-range mash 40, tha later generated by t ha t hro sholdeir 42 aa described 
above, are convolved by convolver 76. The over-range mask 40 is shown if Figa. 4(d) and 9(a) 

The result of this convolution 76 is a trapezoidal signal 88, shown In Fig. 8(b), having half* tha additional 
slops S needed to correct (through multiplication) the data of tha clipped channels 46 to tha slope 58 Identif lad 
by the slope calculator 72. In particular, trapezoidal signal 88 has a peak value of one corresponding to soma 
of the In-rartge channels of the density profile 44 and a slope of half of S. as calculated by the slope calculate? 
72 corresponding to some of the clipped channels of the density profile 44. 

The trapezoidal signal 88 Is next clipped and scaled by two dipper/sealers 104 and 108, In general having 
different clipping and scaling values as will be described, to produce a correction mask 78, used to correct 
the clipped channels 46 for the patient contribution, and a blending mask 108 used to blend the values from 
the headholdar modal 102 with the in-range channels. 

Referring now to Figs. 3 and 9, the correction mask 78 Ea produced by the first dipper/sealer 104 which 
clips the trapezoidal signal 88, shown in Fig. 9(b) produced by the convolver 76 to values less than Pl where 
Pi , ta preferably equal to one-hatf, and than multiplies this clipped signal by 1/p, producing the correction mask 
78, shown in Fig. 9(c), having a valua of one for all of the in-range channels of the density profile 44 and a 
slope equal to the slope S for the over-rang© channels. 

The correction masts 78 Is then multiplied by the density prof Be 44, as indicated by multiplier 80. so as to 
correctthe dels of the density profile 44 for the patient contribution to produce the desired slope 58 more close- 
ly approximating the true density profile 48 in a weighted patient contribution signal 118. 

Th ® b 'f ndin 9 mas{ < 1 08 Is Produced by second clipper/sealer 1 06 rich dips the trapezoidal signal 88, shown 
in F.g. 9(b) produced by the convolver 76 to values less than a predetermined value p* lass than one half 
, i !" 8 " blracte lhe8e va,U8S from Ps to th « '"verted signal and then multiplies this clipped and inverted 
slgnd by 1/p, producing the blending mask 108, shown in Fig. 9(d). The blending mask 108 has a value of 
zero for an of the In-range channels of the density profile 44 and a portion 110, determined by p. of the dipped 
channels M> the portion 110 approximating the area of the dipped channels 46 primarily affected by the con- 
tribution of the patient 12 and thus correrted by the correction mesk 78. The value of p Is determined empir- 
ically The magnitude of the slope of the blending mask 108 Is the seme as that of 8 described above 

The blending mask 1 08 Is used to weight the data from the headhdder model 102, effectively zeroing the 
headholder model for the In-range channels, which implicitly include the actual headholder contribution and 
giving I Increasing weight to the headholdar m<)del 102 for channels less likely to have ettenuation by t to 

2i£ T.? T ?" lon ,or lhe palient contributton v**™ b * the 78 

ST? ^ ^ contribution m **** *'aas aa a result of the reliance on the slope 58 and start- 
2522? II 8 . "5? . ™* imP " Cit,y ° ffeet by the actual ^holder contribution. The weighting of the 
headholder model 102 Is provided by multiplier 112. shown In Fig. 3, which multiplies the blending mask 108 
by the data for the headholder model 102 to produce a weighted headholder comrtbutton ™4 

The weighted headholder contribution 114 and the weighted patient contribution 118 are added together 
at summer 11 6 to produce a corrected projection 80. ^ 

«h. iHSii P "X?! 0 t n8 1 ! 0 formuni P | *9«*'y angle* are then reconstructed according to reconstrudlon math- 

ElU %£T*.Z 0S f *■ * ,h8 ° TO 8UCh a,80r,thm ,8 ™"*» ba <* P^ction^rtch 

lhe Fouriertransform of the projedion 90. It will be recognized that thisfBtering is equivalent to a convolution 

t^STZ^ it8 f """to"*** with thei ™™ FouriertransfornJthef BteVa 
arte latic. and that such convolution effectively spreads errors to the portions SOoverthe entire image. As 
mentioned, image artifacts caused by errors in the projections 90 at the periphery of the patient 12 may create 
image artifacts spreading throughout the image. ^ may create 

block (IT ^ r8COnS,rUCti0^ • lndlcated by process block 82 - *» 'mages are avalabta as indicated by process 

"^/"^ .thattheoperations contained with the dotted line 68 of Fig. 3, representing those operations 

SET^JL^ T 89 re f 0n8,ructor **< are all suitable for pipeline processing. The requlremente for exe- 

♦25 ! T ST P Plp8,,n ° f ° rm i8 Simply that ,ha °P eration be * ^'ng performed on the in- 

putted data of the sequential channels of the projection in their unsorted order without the need to search 

£50L J T mU> th8 Channel d8te that h8S P"**"* ■»«« examined. Because the thresholding of 
block 42 requires only a comparison of each channel of Information againstthe ADC maximum 39. the thresh- 

22 m ^ s , TO y ^ P'Pf "ad.Thesame is truewlththesummation of process block74. and the convolution 
and multiplication and addition of process blocks 76, 80, 112 and 116, Likewise this is true for the clipping and 
scaling of process blocks 104 and 106. a 

The particular Implementation of the correction process employing the calculated slope 58 is thus well 
suited for use in tomographic imaging systems which typically have such processors for performing the recon- 

19 



struct Eon of the images from the projection data. 

Many modifications and variations of the preferred embodiment which will stB be within the spirit and 
scope* the Invention will be apparent to those of ordinary skill in the art For example, the length and shape 
of the box earkornel used in theconvolution may be varied to accountfordlffarentmodels of the Imaged object 
Oearjr tha invention Is not Gmlted to use in correcting ovar-range data associated with beadho.de* but S 
^?s Mrre - Ct f °' ,he attenuatton of anv invari ant within the fan beam. In order to apprised 

puWta of the vanoue embodiments that may faS within the scope of the invention, the following claims are 



Clolomo 
1. 



t °^T * n atCUll ! Ct USS h a ,orn °fl ra Phic scanner, the scanner employing a fixture (15) adjacent 
TJ J"!"T? (12) and mftfl8urin 8 a set of intensity values (36) each related to the Intensity of 

^Tf "f (24) 8,009 ° ne 01 8 8at of adJacent < 21 ' throu Sh both the scanned body and 
the fbdure from a given gantry angle (8). the intensity values indicating the Intensity of the received ra- 
diation up to a maximum (39). and not indicating the Intensity of the received radiation above that maxl- 
mum. the set of intensity values being combined with other sets of intensity values and reconstructed to 
produce an image (86). the circuit comprising: 80 
a threshold detector (42) receiving the set of intensity values for Identifying the set of intensity vat- 

HI! i """T . ten8ity V8lU8S corra8 P ond,n fl to rec ^ed radiation not above the maximum, and over- 
range intensity values corresponding to received radiation above the maximum- 

r«„„ I*" 9 ,* Mm0t (102) reC8Mnfl tha valua of th8 8 antr y engle <°r Producing estimated over- 
[h^i ? correspondin 9 to tha "^sity of radiation along rays associated with those values 

through the fixture alone at the gantry angle; and 

ty^J ° 0mDindr 74 ' n > 7fl - 106 - 112 - ™> ^ combining the in-range intensity values with the es- 
timated over-range intensity values prior to reconstruction of the Intensity values Into an image. 

2. The compensation circuit of claim 1 wherein the threshold detector includes a comparator for comoarino 
m ? f ?\ V£ "V 8 to thB ' naKimUm to ldentlf y only those Intensity values breaterthanjejaualto the max- 
imum Intensity value as over-range intensity values and the rest as in-range intensity values" 

3. The compensation circuit of claim 1 wherein the model generator (102) includes; 

and 8 tobl9 r8C9iVi " S 9Bntry a " fl,e W f0r producin 9 8 Predetermined set of coefficients; 

a function generator receiving the coefficients from the look-up table for generating the estimated 
over-range Intensity values from a polynomial employing the coefficients with e basis fu nrtion (1 00* 

4. The compensation circuit of claim 3 wherein the basis function (100) is piece-wise Gaussian. 

** Llnf.TT" 0 ? ^ ° fdaim3 Wherei " lhe ""^iner replaces a predetermined portion of theover- 

IZ-^Stl ."J tha eStimat9d ° Var - renfle V8,uaa - tha l^ton being those over-rangeva^ues 
associated with rays at the outer edge of the set of adjacent rays. 

8. The compensation circuit of daim 3 wherein the combiner comprises: 

-1 m8Sk ? anefator <f > tor Producing a multipoint binary signal having a first value associated with 
each ,n-range intensity value and a value of zero associated with eech over-range intensity Sue 

ri «r«, Jn T P n l rat ° r {74 ' 72) for Producing a multipoint binary signal having a first value for'a pre- 
determined number of points and a value of zero for all remaining points- ^ 

car 0 J^Z°1 ^J!,^ fW COnVOl ! in8 J 5iflnal from the mask Senerator signal with the signal from the box- 
car generator to produce a convolved signal; 

... nf tt^W? f " mLihip ' yln9 lha eetimated over-range intensity signals by a corresponding val- 
ue of the convolved signal to produce weighted estimated over-range intensity values; and 

a summer (116) for adding the weighted estimated over-range values to the in-range intensity val- 



^Z^SSZ^ ? £ ^ Wharei " ,h8 f ir8t V8,Ue 6qual8 ° nfi and herein the convolver 
2i £ ? flnal ^ th4> m8Sk 8enerator 8i 9 nal afld signal from the box-car gen- 

erator to values less than a predetermined blending value fo) and then multiplies this clipped signal by 
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(he inverse of the predetermined blending value to produce the convolved signal. 

8. In a tomographic scanner measuring a set of Intensity values (36) each related to the intensity of received 
radiation (24) along one of a set of edjacent rays (21, 22) through both an Imaged object (12) and a fixture 
(15), from a given gantry angle (9), the fixture being adjacent to the body to be scanned, the intensify 
values Indicating the Intensity of the received radiation up to a maximum (39). and not indicating the tiv 
tensity of the received radiation above that maximum, the set of intensity values being combined with 
other sets of intensity values and reconstructed to produce an image (86), a method of correcting for over- 
range comprising: 

generating (42) an over-range mask (40) Identifying the set of intensity values as In-range intensify 
values corresponding to received radiation not above the maximum, and over-range intensity values cor- 
responding to received radiation above the maximum; 

determining (102) estimated over-range intensify values corresponding to the intensity of radiation 
along rays associated with those values through the fixture alone at the gantry angle; and 

combining (112, 116) the In-range Intensify values with the estimated over-range Intensity values 
prior to reconstruction of the intensify values into an image. 
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(g) Reduction of Image artifacts from support structures in tomographic imaging. 

@ An apparatus for reducing Image artifacts 
caused by over-ranging or dipping of the data 
collected In a tomographic scan using a patient 
support (15.53), models the density profile of 
the patient support with a polynomial. Coeffi- 
cients of the polynomial are stored in a look-up 
table according to gantry angle 6 and used to 
estimate the over-range data, The estimated 
data is blended with the in-range data by con- 
volving an over-range mask (40) indicating 
which data is in-range and which data is over- 
range with a box car convolution kernel (75) 
which produces a trapezoidal weighting mask 
This weighting mask is multiplied by the esti- 
mated data and then summed to the in-range 
data to correct for the over-range portions of 
the data collected. 
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